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A quantum Hamiltonian is constructed for SU(3) lattice QCD entirely from color triplet 
Fermions — the standard quarks and a new Fermionic "constituent" of the gluon we 
^ call "rishons". The quarks are represented by Dirac spinors on each site and the gauge 

Q , fields by rishon-antirishon bilinears on each link which together with the local gauge 

■ transforms are the generators of an SU(6) algebra. The effective Lagrangian for the path 

' integral lives in i?'* X 5^ Euclidean space with a compact "fifth time" of circumference 

' (/?) and non-Abelian charge (e^) both of which carry dimensions of length. For large 

/3, it is conjectured that continuum QCD is reached and that the dimensionless ratio 
= / (3 becomes the QCD gauge coupling. The quarks are introduced as Kaplan 
chiral Fermions at either end of the finite slab in fifth time. This talk will emphasize the 
gauge and algebraic structure of the rishon or link Fermions and the special properties 
that may lead to fast discrete dynamics for numerical simulations and new theoretical 
insight. 



1 Introduction 

K ' 

0^ ' A quarter of a century after its discovery, solving QCD remains one of the ma- 

jor challenges for theoretical particle physics. While Quantum Chromodynamics 
is generally acknowledged to be a complete theory of all hadronic or nuclear in- 
teractions, only special perturbative consequences are well understood. In 1974 a 
non-perturbative formulation of QCD was given by WilsonEl using a lattice regulator 
. which maps the quantum field theory onto a classical statistical mechanics problem 

in a 4-d Euclidean space. Although a variety of techniques have been developed 
to solve lattice field theories, at present the most powerful tool is the Monte Carlo 
sampling of the partition function of the corresponding classical statistical mechan- 
^ . ics system. Moreover, the most efficient numerical algorithms for QCD suffer from 

critical slowing down when the continuum limit is approached and thus overwhelm 
even the fastest supercomputers. Consequently it is reasonable to look for new 
formulation of the QCD problem. This talk will present a radically new approach 
based on recent work by Browcr, Chandrasekharan and WieseH. Additional back- 
ground material for this approach has already been given by Uwe-Jens Wiese in his 
talktl just preceding this one. 

Here we present our alternative non-perturbative approach to QCD in the 
framework of quantum link models, which leads to a new computational frama- 
work we call the QCD Abacus |[ Such gauge modelsjvere first discussed by Horn □, 
and studied in more detail by Orland and Rohrlich u and by Chandrasekharan and 



"Lecture at "APCTP-ICTP Joint International Conference '97 on Recent Developments in 
Non-perturbative Method" May, 1997, Seoul, Korea. MIT Preprint CTP 2693. 

''Abacus: an instrument for performing calculations by sliding counters along rods. From Greek 
abax (afiai), literally slab. 
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Figure 1: Dimensional reduction in slab geometry: An effective 4-d lattice theory results with 
lattice spacing /3c and a gluonic correlation length l/m which is exponential in (3. 



Wieseu. In these models the 4-d classical statistical mechanics problem of standard 
lattice gauge theory is replaced by a problem of 4-d quantum statistical mechanics. 
In particular, the classical Euclidean action is replaced by a Hamilton operator. As 
a consequence, the ordinary c-numbers in the standard formulation of lattice gauge 
theory for the gauge matrices and their local gauge transformations are replaced 
by non-commuting operators acting in a Hilbert space. Thus one may say that 
the geometry of the gauge manifold is replaced by a non-commutative geometry. 
Nonetheless as discussed in detail in Wiese's talk, both 4-d non-Abelian quantum 
link gauge theories and 2-d quantum sigma models are expected to "dimensionally 
reduce" to an effective continuum field theory exhibiting asymptotically free scaling 
as the extra compact dimension becomes large (see FigJ^). These crucial physical 
motivations vdlllor the most part not be repeated here so the reader is referred to 
the literatureS'Dij for them. 

In spite of these interesting developments, until the recent paper by Brower, 
Chandrasekharan and WieseH, no candidate quantum link theory for QCD itself 
had been proposed due to the difficulty in finding a mechanism to break the U(3) 
gauge group down to SU(3). All earlier constructions suffered from an extra U(l) 
gauge symmetry with a strongly coupled "photon" mode which made any serious 
application to QCD unlikely. Here we describe the algebraic insight that allowed 
us to construct a candidate quantum link Hamiltonian for QCD. The structure of 
our proposed QCD quantum link Hamiltonian is elegantly expressed in terms of 
new Fermionic triplet /anti-triplet fields (3,3) which in a tight binding mode form 
"composite" gluonic degrees of freedom. We call these Fermionic constituents of the 
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gauge links, "rishons" or link FcrmionsLO. We believe this quantum link formulation 
of non-pcrturbative QCD is both correct and compelling for its simplicity. Perhaps 
it will lead to more powerful methods for attacking long-standing problems in lattice 
gauge theories by analytical and/or numerical methods. 

This talk is organized as follows. In section 2 quantum link models are formu- 
lated in terms of rishons — the Fermionic constituents of the gluons. The inclusion 
of quarks is discussed in section 3 within the framework of Shamir's variant of Ka- 
plan's Fermion method. Section 4 contains a discussion of the prospects for deeper 
theoretical understanding, an analytical approach to the large N limit, for example, 
and cluster methods for Monte Carlo simulations of quantum link gauge theories. 

2 The QCD Abacus 

The main goal of this talk is to display the rules for the quantum link Hamiltonian in 
a form that allows one to appreciate the algebraic structure in a vivid and intuitively 
appealing form. As illustrated in Fig. ||, the gauge sector of the quantum link QCD 
Hamiltonian acts on a discrete Fermionic Hilbert space, represented by colored dots 
for occupied Fermionic states associated with each link. The link occupation level 
is set to half filling for the sum of left triplet (3) and right anti-triplet (3) link 
Fermions. The plaquette term in the Hamilton operator slides the beads in the 
Abacus along the link in the clockwise or anti-clockwise direction changing the flux 
by one quantum. Gauge invariance demands that the same color that leaves a site 
must move into that site from the adjacent link in the plaquette, conserving rishon 
number for each color at the sites. Another term is introduced into the Hamiltonian 
to remove the unwanted U(l) gauge invariance which creates and destroys rishon 
"baryons" at opposite ends of each link, reducing the gauge group from U(3) to 
SU(3). Once the origin of these rules is explained constructing other quantum link 
Hamiltonians becomes a natural game of inventing appropriate "cellular automata" 
rules on the Quantum Abacus. 

At first it may seem hopeless or even paradoxical that one might represent 
exactly a continuous gauge group (or even a global Lie group) with a finite set of 
bits per lattice site. However, as discussed in more detail in Wiese's talk, a quantum 
spin is just such an object. For example the electron at rest has only 2 states but still 
respects exact rotational invariance via quantum superposition. Consequently the 
trick is to consider the appropriate quantum operators in a discrete Hilbert space. 
Just as in all quantum systems, the eigenvalues for sets of commuting operators 
enumerate the states and the operators for Lie group symmetries "rotate" them 
continuously into linear combinations of these eigenstates. 

2.1 U(l) Warm Up Exercise 

Before considering the full problem of non-Abelian gauge theory, it is instructive to 
begin with the simpler compact U(l) Abelian gauge theory. With hindsight many 
of the basic concepts needed for full QCD can be introduced. Since we want to 
construct a 4-d gauge Hamiltonian, it is useful (again with the benefit of hindsight) 
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Figure 2: Rishon dynamics: The trace part of the Hamiltonian induces a hopping of rishons of 
various colors around a plaquette. The determinant part shifts a color-neutral rishon baryon from 
one end of a link to the other. 



to consider a Kogut-Suskind Hamiltonian ^ in 4-d as an example, 

^ - 1 . . . . 

Hks = ^^Ex.4iEx^^,+ ^ ■^Ux,,_iUx+^i,i^Ux+u,J-^l,i^ ■ (1) 

The sums extend over lattice sites x and directions = 1,2,3,4. Local gauge 
invariance on each link (x, x + jj) requires, 

[E,U]^U [E,U^] = -U^ [f/,?7^] = 0, (2) 

where to simplify the notation the x, /i subscripts have been dropped '0. 

It is easy to prove (just as with the canonical commutation relations [q^p] = i) 
that the only solutions to this algebra are infinite dimensional "matrices" . For 
example, in a basis in which all the ?7's are diagonal (the "field basis"), E is given 
by differentiation, 

U^e'^ , ^--4> (3) 



'^If we had considered the Kogut-Suskind Hamiltonian for a non-Abelian theory, we would have 
had to distinguish between the electric fields that act as right generators, B^, and those that act 
as left generators, Ei^ = 0^ Efi(j^ . 

^ Including the link subscripts would lead to equations such as [-Ej, jj, l/y^^,] = Ux,ii.&x,y&ti.,v, 
throughout. Since all commutators for operators on different links are zero dropping these indices 
should not cause confusion. 
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where 9 = eA. On the other hand, if we aUow the commutator [[/, [/t] to be 
non-zero, finite dimensional representations can be found. This amounts to a new 
quantum prescription with the group manifold represented as a non-commutative 
geometry. For example in the smallest representation, the rest of the algebra is 
easily satisfied by 2 by 2 sigma matrices with 
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and the gauge generator, 



E 



5^3 



-1 



(4) 



(5) 



(6) 



In a basis where the E field is diagonal (the "electric basis"), we may view this 
2 by 2 representation as simply a truncation of the infinite dimensional represen- 
tation to a single unit of fiux. To show this, let us introduce a set of Fermionic 
creation/destruction operators. 



{c;,4} = 



one for each Fourier mode I = ±i, ±|, ±|, • • • , ±cx), and constrain the Hilbert space 



to 



J2 
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Now the U{1) commutation relations (g) are satisfied by the operators, 

U - / ^0(0)e*^0t(0)^^c,4^i, (7) 

(8) 



where 



AW 



It is now straight forward to demonstrate that the sole 



consequence for the algebra of truncating the fiux ( i.e. restricting the sum in 

L-l 



= aie^^^ to |Z| < -/^ ) is to introduce the non- vanishing commutation relation. 



It is important to contrast this truncation procedure in the "electric basis" relative 
to the more conventional approach of replacing the group manifold by a finite set 
of angles 6'„ — nA9. In the latter, the continuous symmetry group (i.e. U{1)) is, 
at best, replaced by an approximating subgroup (i.e. Z^), whereas by truncating 
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to low dimensional representations in the Fourier space the full symmetry group is 
left intact! 

By taking the lowest non-trivial flux on a link {I ~ i^), we obtain the smallest 
representation, 

U = afet [/t = _at5 R = b^b a^a (9) 

where a = c_i,b = ci. These satisfy the commutators 

[R, U]^U , [L, U]^^U , [U, U^=R-L . (10) 

In the 2 dimensional subspace with a'^a+b'^b = 1 and E = {R~L)/2, this reproduces 
exactly our earlier sigma matrix construction in Eqs. (Q^). 

Let us summarize the situation. We can construct a U(l) gauge invariant 
Hamiltonian acting on a 4-d lattice simply by replacing the link phase factors, 
exp[±i0], by sigma matrices, . The electric field term, being the value of the 
quadratic Casimir for the lowest dimensional representation, is just an additive 
constant and can be dropped. The resulting quantum link Hamiltonian for U{1) is 

^ = '^i,i^'^x+t,,^(^x+u,^,(^x,v (11) 
and its partition function is 

Z = Trexp(-/3F). (12) 

We note that in this Abelian example, we have made use of two equivalent repre- 
sentation: one in terms of explicit (sigma) matrices and another in terms of link 
Fermions (or rishons). The strategy in the next section is to explore both of these 
approaches to generalize this construction to SU(N) Yang-Mills theory. 

2.2 U(3) Quantum Link Theory 

Let us consider a quantum link Hamiltonian, 

'^^\K^^U^+f^,^lJl^i>.^PU + ••• (13) 

for the non- Abelian field operators Ul? ^- In this Hamiltonian each plaquette is the 
trace of the matrix product in color indices i, j — 1, ■ ■ ■ , N , and a direct product of 
q number matrix elements. The dots in Eq. |l^ represent other operators that will 
be needed for full QCD, for improved actions, etc. 

For SU(3) the local gauge algebra is realized as an SUr{3) x SUl{3) symmetry 
for each link operator U, 

[RaiUtj] ^ UikXkj , [La,Uij] ^ -Xfi^Ukj , (14) 

where for notational simplicity we again suppress the link indices x, ^ in the above 
formula: U = Ux.fj., R = Rx,tj. and L = L^^fj,- Note that each link in fact has 9 oper- 
ators Uij for I, j = 1, 2, 3 transforming as a (3, 3) representation Ij. All commutators 

■^A better notation for these tensors would be to use upper and lower indices C/^ but the standard 
practice in lattice gauge theory ignores this nicety. 
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involving different links are zero. The SU(3) generators obey the Lie algebra, 

[Ra, Rp] = 2ifafjyRj , [La, Lf}] — 2ifa/3^L^ , [i?Q,L^]=0, (f5) 

where the Gell-Mann matrices ■[] satisfy the algebra, 

AaA/3 ~ 5al3 + ifa/Si^j + daf3jXj ■ (16) 

If we introduce infinitesimal generators, 

for local SU{3) gauge transformations, it is trivial to demonstrate that as a con- 
sequence of the above Lie algebra the single plaquette Hamiltonian j ]!^ ) does obey 
local gauge invariance. Thus we have a general algebraic framework for construct- 
ing quantum link Hamiltonians with local gauge invariance. Surely this framework 
is not unique, but given the reasonable assumption that all non-zero commutators 
are confined to operators on a single link, it is probably the most natural. 



2.3 Minimal Solution 

As we will explain in more detail in the next section, there are many solutions 
to this algebraic problem. Here we attempt to construct the lowest dimensional 
non-trivial solution, which for U(N) gauge theory will turn out to be built on the 
smallest irreducible representation of SU(2N). This representation has dimension 
2N. 

We introduce the ansatz. 



1 





-iAa 




(18) 



and 



Ro 







Lr 





A, 



(19) 



It turns out that the gauge algebra (|lj) is satisfied, if the q-number link operators 
Ac obey exactly the same Gell-Mann algebra as the c- number generators Aq, 



AqA^ — 5al3 + ifalSjAj + da/SjAj . 



(20) 



Thus the A's may also be represented as 3 by 3 matrices for each link. To prove 
this, substitute the ansatz into Eqs. (Ill) to get the two relations. 



Rail + XpA,}) 
{l + X0Ap)La 



(1 -I- Xi}Af})Xa 
Xa{l + XpAp). 



(21) 
(22) 



^To avoid annoying factors of ^, I have absorbed this factor into redefining the A 's and the 
f and d structure constants. This corresponds nicely with the traditional convention for SU(2): 



Oq/3 ■ 
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These are satisfied identically. For example expanding the first relations (21) above 
yields 



Aa + + i/Q^-yA/3A-y + dct^^A^A^ — Aq + Aq + i/^Q-yA^A^ + d^Q^A/jA^ . (23) 

Using only the condition that fa/i-y = Tr[[Xa, Xp], X^]/6i and dafi^y = Tr[{Xa, A^}, A^]/6 
are totally anti-symmetric and symmetric tensors respectively, the identity holds. 
The second relation ( p2|) is the conjugate of the first. (Clearly the same construction 
will work for any Lie group.) 

In addition to the SU(3) gauge transforms discussed above, there is a pair of 
left and right U(l) gauge generators. 



Rn — 



1 



















, Lq — 





1 



(24) 



The sum, J\f — Rq + Lq, is a trivial symmetry since it commutes with the link 
operators and the other generators, but the difference, T = Rq — Lq, leads to an 
additional (unwanted) gauge symmetry, 

[T, - 2(7,,^ . (25) 

The consequence of this extra [/(I) generator is that the standard plaquette interac- 
tion in our Hamilton operator (|l^) actually represents a U (3) gauge invariant yaog 
Mills theory. Indeed because of this unwanted symmetry, all earlier attempts alj'Ll 
did not succeed in constructing a Yang-Mills Hamiltonian suitable for QCD with 
only the SU(3) gauge group. 

One way to understand the problem is first to consider the standard single 
plaquette {Tr[Up]) Wilson action and to note that this expression is also invariant 
with respect to U(3) gauge transformation. Consequently, the restriction to SU(3) 
for the Wilson theory is in fact introduced "by hand" by defining the path integral 
with the Haar measure dU restricted to the SU(3) manifold. Therefore, in the 
quantum case the analogous procedure is to insert a symmetry breaking operator 
into the trace similarly restricting the phase space. A natural approach (which 
we will ultimately implement) would be to add to the Hamilton operator a term, 
J'Det[Ux,fj,], for each link. But in our present minimal construction, the direct 
product of a link with itself vanishes ( W^U''^ = 0), so such a determinant term 
vanishes identical as an operator! The problem can be traced to the fact that the 
smallest representation just does not allow one to squeeze more than one unit of flux 
onto a single link. One (rather inelegant ) way out is to build the SU(3) determinant 
out of three different operators spreading the flux over three different paths joining 
the point x to x + fi. Fortunately, the rishon representation will afford a more 
elegant and minimal solution to constructing the determinant locally on each link. 



2.4 The Rishon Representation 

We now proceed to this approach to constructing quantum link QCD based on 
Fermionic triplets (or rishons), which naturally leads to higher dimensional repre- 
sentations of SU(6). The basic algebraic structure of these rishon models is exactly 
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the same as our minimal model given above. In particular all the commutators (and 
therefore the group properties) are unchanged. 

In addition to the gauge commutators (^4|), let's enumerate the other com- 
mutators. Again the only non-zero commutators for Ux,fi and ^ involve the 
same link {x,x + ^) so we suppress the link index in discussing the algebra. With 



U = Uo + iKUo, and C/t = ul - iXa&l, we have 

[Uo,U^o] = T, (26) 

[Ua^U'p] — SapT + ifap^{Ry + L^) + dai3^{R~f — L^) (27) 

[Uo,Ul] = l(i?a~£a). (28) 

Also the commutators, 

[U^,Up]^0 , [C/t,J7t] ^0, (29) 



are zero for all a, (3 — 0, ■ ■ ■ N"^ — 1. In addition in the minimal construction above, 
we also have the nil potency property that lJx,^^ and ^ just as in the U(l) gauge 
theory act like a creation-destruction operators: Ux,ii = ^^^^ — 0. Not being part 
of the Lie algebra, we will see that this special nil potency condition is a peculiarity 
of the lowest dimensional representation. 

The rishon construction begins with the following ansatz. Consider 2A^ Fermionic 
operators ai and hi for i = 1, ■ ■ ■ , N , with anti-commutators, 

{ai,a]}=5i.j , {bi,b]}^8ij, (30) 

on each link and all other anti-commutators zero between different links and be- 
tween a's and b's on the same link. The local gauge rotations are 

^a-E«I^S«^- ' Rc.^Y.^\\%b, , (31) 

i i 

on each link with a = 0, 1 • ■ • N'^ — 1. Just as in traditional current algebras, these 
bilinears obviously satisfies the correct SUr{N) x SUl{N) gauge algebra. The 
gauge fields are then represented by 

Uij = E A"j anXrim^lnJ (32) 

or 

U^J = a^b] , (33) 

since the N'^ A's form a complete hermitian basis. Again it is straight forward to 
see that the complete set of bilinears made from 2A^ Fermi fields must give the 
Lie Algebra relations for SU{2N). There is a very appealing interpretation. The 
bosonic gauge field is a Fermion-anti-Fermion bound state. Like the tight binding 
Hamiltonians of solid state physics, these "rishons" constituents never move relative 
to each other, so there are no new internal degrees of freedom. 
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The rishon bilinears act on a 2^^ dimensional Fock space at each hnk which 

can be enumerated by occupation numbers 0,1 for each Fcrmionic mode: aj,6j. 
Moreover the entire algebra commutes with the link number operator, 

7\r=^o + ^o = ^alai + 5^6l6i , (34) 

i i 

SO that the 2^-'^ dimensional matrices are reducible into a set of 2A'' + 1 irreducible 
representations. The enumeration of the representations go as follows. The link 
number Rq + Lq = K extends over K — 0,1, ■ ■ ■ 27V, but there is a Z2 particle- hole 
symmetry that relates K to 2N — K so really the models are distinguished only 
for K = 0,1,- --N. The case K = is trivial. Consequently, there are only N 
distinct representations. For U(l) there is a unique model which is exactly the one 
discussed above in Section 2.1. For general N, the case K = 1 corresponds to a 2N 
dimensional representation, since the Fock space is restricted to a single Fermion 
with one of the a^'s or one of the fe^'s occupied. This is exactly the minimal U(N) 
construction presented above. The models with K > 1 are new. For SU(2), K = 2 
is a 6 dimensional representation of SU(4) (anti-symmetric two index tensors). The 
square of a link field now has one nonzero component given by 

det[U.,^.] = l^ eij ei^yUl;:^^^;, = 2! aH^Jjljl^. (35) 

Consequently a term which breaks U(2) down to SU(2) can be introduce locally on a 
single link in the 6 dimensional representation. For SU(^3) there are two new models 
— for ii' = 2, a 15 dimensional representations with U'^ non-zero but C/^ = and 
for = 2. a 20 dimensional representations with C/^ non-zero. Just as in SU(2), 
the 20 dimensional representation allows a U(l) breaking term to be introduced on 
a single link, 

detp.A = ^ ^^jk U'rw tJ^filiit}^ = 3! aliyj^yjjljljl^ . (36) 

The intermediate case for SU(3) with K = 2 will allow U(l) breaking to be intro- 
duced on a single plaquette coupling a di-rishon with a rishon. Clearly there is a 
tradeoff between the size of the representation and the ability to represent operators 
locally on the lattice. 

In summary, we pick our candidate Hamilton operator for SU(N) Yang-Mills 
theory to be 

x,n 

The plaquette part of the Hamilton operator shifts single rishons from one end of a 
link to the other and simultaneously changes their color. The determinant part, on 
the other hand, shifts an entire "rishon-baryon" — a color neutral combination of 
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rishons — along the link. Each link is restricted to a half filled state of N rishons 
by fixing, 

so that the link algebra carries a 2A^!/(iV!A^!) dimensional irreducible representation. 
This rather simple rishon dynamics may facilitate new analytic approaches to QCD, 
and may also be useful in numerical simulations of quantum link models as discussed 
in Section 3. 



2.5 Comments on SU(2) Gauge Actions 

SU(2) waSpthe first example of SU(N) Yang Mills theory formulated as a quantum 
link modela. For SU(2) there is a special simplicity, which allowed the link operators 
tj^ to be represented by 7^ matrices in an S0(5) algebra, but this special feature 
tended to obscured further generalizations. With our general rishon method for 
SU(N), it is natural to ask if this special S0(5) construction for the SU(2) model 
represents a distinct alternative. In fact this is not the case. Following the rishon 
construction, one can develop the SO (5) algebra for SU(2) as a limiting case. 

Returning to our minimal construction presented in Sec. 2.3, the algebraic 
distinction between SU(2) and SU(N) is due to the d-symbol. To remove the extra 
U(l) gauge invariance for SU(2), one makes use of another form of the link operator, 
Vij = {a2U^(J2)ji that also obeys the same gauge algebra for SU(2), where U is given 
by the expression in Eq. This is possible for SU(2) because there is an extra 
G parity symmetry: CTq ^ 0'2<Ja<^2 = ^ = — that preserves the algebra. 
(Note invariance of aaap = Sa(3 + i€af3-y<^'y under CTq — (cq)"^.) In other words for 
SU(2), the (2) and (2) representations, which are related by the raising and lowering 
symbols ( = e*-' ~ 'i{<y2)ij), are unitarily equivalent. 

However if we attempt the same construction for SU(3) by introducing. 



V = 

















1 





iAa 






(38) 



the new fields fail to obey the gauge algebra (|lj 

-{l~XpA0)Ro, ^ il-XpAp)Xc, , (39) 
La{l~Xf3Ap) 7^ -A„(1-A/3A^) , (40) 

due to the opposite sign for the d-symbols on the right and left hand sides. Conse- 
quently only for SU(2), do we have the special construction, 

Hsu(2) = ^ E Tr[W,,^Z^.+^..Wt^,_^Z^t J , (41) 

x;p.>iy 

where the new link operators are a sum of two terms, I4x,f_i = Ux,ij. + Note also 
that U =U{1 + T) /2 'AndV = U(l- T) /2. This Hamihonian is "by accident" both 
hermitian and gauge invariant under SU(2) but not U(2). 

However another construction is also possible for SU(3), beginning with the 
minimal 6 dimensional representation in Eq. This is based on the observation 
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that any link U{x, y) connecting x to y, transforms as a (3, 3) representation of 
Ur{2>) X J7l(3). Thus not only can we form U(3) gauge singlets by contractions 
with (5ij's to make "mesons" at the sites but we can also form SU(3) singlets by 
contractions with eyfe's to make "baryons" at the sites. For example we can add a 
term, 

Hi^Ji e,jkU^^Ax,y)U^'}{x,y)ui%{x,y)e.,,,k' (42) 

<x,y> 

where the U's are 3 independent paths between x and y. 

In the special case of U(2), one may use two paths connecting x and y on a 
single plaquette. With appropriate weights and paths, this procedure can lead you 
back to the earlier S0(5) suggestions for an SU(2) Hamiltonians (|4l|). For example, 
take 

Hi^Ji E y)^£' y) '^'r ■ (43) 

<x,y> 

Making use of the raising and lowering symbol eij = {i(J2)iji this expression can be 
re-written as 

Hi = Ji E uS^{^,v)vif{y^^) (44) 

<x,y> 

Identifying [7*^^^ and F*^^^ with products of links on the edges of the plaquette, we 
get the individual terms in Eq. ^ that break U(l) gauge symmetry. The former 
choice (^ ) is just one possible Hamiltonian, where the epsilon symbol is inserted 
or not into all the corners of the plaquette with equal weight. 



3 Quark Fields and Full QCD 

To represent full QCD, one must add quarks to the quantum link model. This is 
more or less straightforward- although some subtleties arise related to the dimen- 
sional reduction of Fermionsa. 

Again the guiding principle is to replace the classical action of the standard 
formulation by a Hamilton operator which describes the evolution of the system in 
a fifth Euclidean direction. For the quarks this implies replacing the Grassmann 
variables ipx, i^x by 'I' J, 75, where ^t]. and '^x are quark creation and annihilation 
operators with canonical anti-commutation relations 

,j,,6/3t} ^ 5^y5,,5ab5c.p, {^T . "^ f) = {*r^ - 0, (45) 

where (i,j), {a,b) and (a,/3) are color, flavor and Dirac indices, respectively. 
The full QCD quantum link Hamiltonian is now given by 

+ ^ E [2*h5*. - ^h,Ux,^.^x+^ - *I+/,75?7.\^*J- (46) 
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The generators of SU{N) gauge transformations must include quark bilinears, 

G^ = E(^"-A,. + ^".M) + *i^"*- (47) 

but it is still trivial to show that H commutes with all G". 

To ensure the proper dimensional reduction of the quarks, their boundary condi- 
tions in the fifth direction must be chosen appropriately. The standard antiperiodic 
boundary conditions, which are dictated by thermodynamics in the Euclidean time 
direction, would lead to Matsubara modes, = 27r(n5 + However for our 

5-th time, this would lead to a mass of 0(/3) for the dimensionally reduced quark. 
On the other hand, the confinement physics of the induced 4-d gluon theory takes 
place at a correlation length which grows exponentially with /3. Therefore, quarks 
with antiperiodic boundary conditions in the fifth direction would remain at the 
cut-off freezing out the quarks and returning the dimensionally reduced theory to a 
quarkless pure Yang-Mills theory. One possible solution is obvious. Simply choose 
periodic boundary conditions for the quarks in the fifth direction. This gives rise 
to a Matsubara mode, ps = 0, that survives dimensional reduction. Since the ex- 
tent of the fifth direction has nothing to do with the inverse temperature (which is 
the extent of the Euclidean time direction), one could indeed choose the boundary 
condition in this way. 

There is another perhaps more attractive possibility. The above scenario with 
periodic boundary conditions for the quarks still suffers from the same fine tuning 
problem as the original Wilson Fermion method. The bare mass matrix A4 would 
have to be adjusted very carefully in order to reach the chiral limit. In practice this 
is a serious problem for numerkial simulations. This problem has been solved very 
elegantly in Shamir's variant eI of Kaplan's Fermion proposal Kaplan studied 
the physics of a 5-d system of Fermions, which is always vector-like, coupled to a 
4-d domain wall that manifests itself as a topological defect. The key observation 
is that under these conditions a zero mode of the 5-d Dirac operator appears as a 
bound state localized on the domain wall. From the point of view of the 4-d domain 
wall, the zero mode represents a massless chiral Fermion. The original idea was to 
construct lattice actions for chiral gauge theories in this way. 

Subsequently, Shamir pointed out that the same mechanism can also be used 
to solve the lattice fine tuning problem of the bare Fermion mass in vector-like 
theories like QCD. He also suggested a variant of Kaplan's method that has several 
technical advantages, and that turns out to fit very naturally with the construction 
of quantum link QCD. In quantum link models, we already have a fifth direction for 
reasons unrelated to the chiral symmetry of Fermions. We will now follow Shamir's 
proposal, and use the fifth direction to solve the fine tuning problem that we would 
have with periodic boundary conditions for the quarks. The essential technical 
simplification compared to Kaplan's original proposal is that one now works with 
a 5-d slab of finite size /3 with open boundary conditions for the Fermions at the 
two sides. This geometry limits one to vector-like theories, because now there are 
two zero modes — one at each boundary — which correspond to one left and one 
right-handed Fermion in four dimensions. This approach fits naturally with our 
construction of quantum link QCD. In particular, the evolution of the system in 
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the fifth direction is still governed by the Hamilton operator of Eq. The only 
difference relative to Wilson's Fermion method is that now r < 0. 

The partition function of the theory with open boundary conditions for the 
quarks and with periodic boundary conditions for the gluons is given by 

Z = Tr(0|exp(-/3i/)|0). (48) 

Here the trace extends only over the gluonic Hilbert space of the quantum link 
model, thus implementing periodic boundary conditions for the gluons. The bound- 
ary conditions for the Fermions are realized by taking the expectation value of 
exp{—(3H) in the Fock state |0), which is annihilated by all right-handed '^rx and 
by all left-handed '^]^^ V3. As a result, there are no left-handed quarks at the bound- 
ary at — 0, and there are no right-handed quarks at the boundary at — (3. Of 
course, unlike periodic or antiperiodic boundary conditions, open boundary condi- 
tions for the Fermions break translation invariance in the fifth direction. Through 
the interaction between quarks and gluons, this breaking also affects the gluonic 
sector. We don't expect this to be problematic, because we are only interested in 
the 4-d physics after dimensional reduction. The only crucial feature is that both 
quarks and gluons have zero modes that survive dimensional reduction. 

There is one technical detail. We have argued before that the gluonic correlation 
length is not truly infinite as long as [3 is finite, but — due to confinement — it 
is exponentially large. The same is true for the quarks, but for a different reason. 
In fact, already free quarks pick up an exponentially small mass due to tunneling 
between the two boundaries, which mixes left-handed and right-handed states, and 
thus breaks chiral symmetry explicitly. 

The important observation is that the bare mass parameter M (considering a 
single flavor for the moment) is not the physical quark mass. The physical mass 
of the quark goes exponentially fast to zero in [3 in the continuum limit without 
fine tuning. At the same time, the doubler modes remain at the cut-off. Let us 
first consider the doubler Fermions, which are diaracterized by n = 1, 2, 3, 4 for the 
corners of the Brillouin zone. It can be shownH that the doubler masses diverge, if 
we choose 

M -f 2nr < 0. (49) 

In fact this precise inequality applies for non-interacting quarks with zero gauge 
fields. But close to the continuum limit, we expect nearly the same inequalities. 
Thus there is a region with r < —M/2, for a sufficiently strong Wilson-term with 
an unconventional sign, where the doubler Fermions are removed from the physical 
spectrum. On the other hand, the mass of the physical Fermion (the n = state) 
is 

fj, = \E{6)\ = 2Mexp(-M/3), (50) 

which is exponentially small in /3, again up to small corrections when we turn on 
the gauge fields in weak coupling. These results indicate how one may avoid the 
fine-tuning problem of the Fermion. The confinement physics of quantum link QCD 
in the chiral limit takes place at a length scale, 

^°^^"P((ll^-27V,)Are2)' (51) 
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Figure 3: Dimensional reduction in slab geometry: The gluonic correlation length 1/m is exponen- 
tial in /3. Left and right-handed quarks live on the two sides of the slab. The overlap of their wave 
functions n and ^I/j^ induces a mass fj,. The corresponding correlation length l/p, is exponentially 
larger than 1/m. 



which is determined by the 1-loop coefRcient of the /^-function of QCD with Nf 
massless quarks and by the 5-d gauge couphng e. As long as one chooses 

247r2 

^ > iUN-2N,)Ne^ ^ ^''^ 

the fictitious induced mass for the chiral quark is exponentially smaller than the 
QCD scale Aqcd and there is a window in which the chiral limit is reached au- 
tomatically as one approaches the continuum limit taking /3 large. For a given 
value of r one is limited by M < —2r (note that r < 0). On the other hand, one 
can always choose J (and thus e^) such that the above inequality is satisfied. The 
corresponding geometry is shown in Fig. ^. 

Of course, we also want to be able to work at non-zero quark masses. Following 
Shamir, we do this by modifying the boundary conditions for the quarks in the fifth 
direction. Instead of using 2:4, 0) — ^Jnix, x^, (3) = we now demand 

Tfl Tfi 

iPl{x, X4„0) ^ ~—^Pl{x,X4,P), tPr{x,X4,I3) ^ -—-ipRix.XA^O), (53) 

where m is a mass parameter. This reduces to the previous boundary condition 
for TO = 0, while it corresponds to antiperiodic boundary conditions for m = 2M, 
and to periodic boundary conditions for m = — 2M. Solving the above equations 
with the new boundary condition indeed yields physical quarks of mass to in the 
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continuum limit /? ^ cx3 as long as m ^ M . It has been shown in Ref. ^ that in 
the interacting theory m is only multiplicativcly renormalizcd. On the level of the 
partition fimction, 

Z = Tr[exp(-/3i7)0(m)] , (54) 
the new boundary condition manifests itself by a mass-dependent operator, 

0{m) = n(*«-*k + *L^*i?.)(*L*L. + *L.^*L), (55) 

X 

which was constructed by Furman and Shamir 0. Note that in Eq. ^ the trace 
is both over the gluonic and over the Fermionic Hilbert space. In the chiral limit 
m = 0, the operator 0{m) reduces to a projection operator |0)(0| on the Fock state 
introduced before. For m = 2M , i.e. for antiperiodic boundary conditions, the 
operator ©(rn) becomes the unit operator, and the partition function reduces to 
the well-known expression from thermodynamics. For m = —2M, i.e. for periodic 
boundary conditions, the operator 0{ra) becomes an alternating sign (—1)^ where 
F is the total fermion number and the trace becomes a super trace. 



4 Concluding Remarks 

In this brief talk, I have attempted to outline the algebraic steps required to convert 
a standard lattice gauge theory into a quantum link theory for the QCD Abacus. 
This has resulted in a demonstration that we can in fact construct a 4-d quantum 
link Hamiltonian for quarks and gluons which obey all the constraints of local gauge 
invariance and discrete symmetries for QCD. There has not been time in this talk for 
many other important theoretical issues such as why we believe this is a legitimate 
alternative formulation of lattice QCD or time to report on the substantial progress 
being made to convert our new scheme to a viable theoretical and computational 
tool. In a sense for quantum link QCD, we now have to repeat the last 20 years of 
work that has been done in the Wilson formulation. 

Although it has not been emphasized, it should be plausible that our general 
approach to converting a classical action into a corresponding Hamiltonian operat- 
ing on a discrete space in one extra dimension can be applied to a wide range of 
quantum field theories. In fact we have coined the term "D-theory" for the gen- 
eral class of quantum spin and quantum link theories JJaat are obtained by such a 
discrete "re-quantization" of path integral expressions Eil. In some ways the proce- 
dure in D-theory is analogous to the construction of M-theory whereby the string 
co-ordinates for the world sheet embedded in 10 dimensions are replace by matrices 
living in 10-1-1 dimensions. When the extra dimension is compactified, the ratio 
(e^//3) of the 11-d coupling (e^) and the compact circumference (/3) becomes the 
coupling {g^) in the 10-d dimensionally reduce theory: e^//3 g^. The general 
procedure in D-theory is to study the symmetry properties of the target field the- 
ory and to replace the co-ordinates (i.e. fields) with a non-commuting algebra that 
expresses all the basic symmetries. Once the representations for the fundamental 
quantum spins (matter fields) and quantum links (gauge fields) have been chosen, 
a large variety of Hamiltonians can be defined. If there is a zero gap in the d-|-l 
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dimensional D-theory before compactification of the extra dimension, then we be- 
Heve that dimensional reduction to and universality with the target field theory is 
generally assured. 

The use of anti-commuting fields to build finite dimensional representations of 
symmetry groups (i.e. rishons in our terminology) is of course a widely used de- 
vice. Indeed the rishon formulation resembles th£_^Schwinger boson and constraint 
Fermion constructions of quantum spin systems E£l. Just as it has for spins models 
this trick may provide new theoretical insight leading perhaps to new analytic ap- 
proaches to QCD, for example, in the large N limit. In particular we note that the 
Hamilton operator (|6|) of quantum link QCD with a U (N) gauge group (i.e. drop- 
ping the determinant term) can also be expressed in terms of color singlet glueball, 
meson and constituent quark operators. These objects consist of two rishons, two 
quarks and a quark-rishon pair, respectively. To see this, just contract the color 
triplet indices in Eq. defining color singlet bilinears for the new fields, 

^x,+pi.+u — C^.fj.C^.u 1 ^^^x — ^x ' '^x,+^l — ^x.ji^x J 1^0; 

representing glueball, meson and constituent quark mesons respectively. In this 
formula, we have defined c^_^ = ax,^^ and c^,_^ = 6^^^ for = 1, 2, 3, 4. The 8-|-4A^/ 
triplet fermion operators obviously combine as bilinears to give the generators for 
an SU{16 + 8Nf) algebra. In principle this theory can be bosonized in terms of 
these degrees of freedom. If large N techniques can be applied successfully, major 
progress in the non-perturbative solution of QCD and other interesting field theories 
is to be expected. In any case, it is intriguing that U{N) quantum link QCD has two 
representations — one exclusively in terms of colored quark and rishon Fermions, 
the other in terms of color-singlet bosons. Returning to SU(N) only adds rishon 
baryons to the set of color singlet fields. 

Although it has not been emphasized in this talk, one can argue ErHu that there 
is a reasonable scenario for this theory to undergo dimensional reduction so that 
the continuum limit is reached for large values of the 5-th time extensions (3 and 
that the resulting theory is universally equivalent to Wilson's formulation of QCD. 
Let us remind ourselves of the analogous situation for the conventional Wilson for- 
mulation. Here the "demonstration" that lattice QCD gives the correct continuum 
theory involved re-deriving renormalized perturbation theory from the lattice and 
showing numerically that the confined phase extends to zero (bare) coupling with 
asymptotically free scaling. In spite of numerous technical difficulties, there is a 
general consensus that these tests have been met for the Wilson approach. While 
we have not yet accomplished a similar degree of certainty for the quantum link for- 
mulation, we believe it can be accomplished. For Quantum link QCD, perturbation 
theory is more difficult, but we have begun to consider a coherent state formalism 
which in principle will allow one to define "spin waves" (i.e. gluons) around a weak 
coupling vacuum. These spin waves should interact in weak coupling like ordinary 
QCD perturbation theory. 

One approach is to begin with the equations of motion for the quantum link 
Hamiltonian, 

4:Ux,^L = i[H,IJx,ii] , ^fx,^ = i[H,fx,f,] 
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= ^[H,Rx,p] , -^L'^^^ - i[H,L'^^^] . (57) 

and argue that in the appropriate classical limit, they lead as expected to the 5-d 
Yang-Mills equation for the mean fluctuations , 

A^{x) = ^ < U{x + fi,x) - U^x + fi,x) > , (58) 

about the vacuum < U >= 1 in the "5-th time" ^5=0 gauge. In addition, 
we are investigating the ph£pi_atructure of 5-d QCD to establish the existence 
of a gapless Coulomb phase EjliJ. The mechanism we postulate for dimensional 
reduction of gauge theories is that compactifying the 5-th_time leads to a gap and 
confinement. As argued by Wiese and Chandrasekharan El, this follows closely an 
analogy with dimensional reduction of the quantum spin Hamiltonians for 2-d non- 
linear sigma models, where there is also a target asymptotically free theory and the 
gap is guaranteed by the Mermin-Wagner-Coleman theoremtJ. Interestingly, in a 
non-Abelian gauge theory, the analog of the Mermin-Wagner-Coleman theorem is 
confinement itself. The direct demonstration of this mechanism for quantum link 
gauge theories (since no one has succeeded to date in defining the confining vacuum) 
will probably come from numerical simulations, where in the slab geometry one 
identifies asymptotic freedom emerging as P is increased. This will require efficient 
Monte Carlo algorithms. 

There is one last issue of tremendous promise and importance that must be 
mentioned. In simulations, the potential advantage of the quantum link approach 
is the application of cluster algorithms as typified by the classic Swendsen-Wang 
algorithmll3 that hfis revolutionized simulations for the Ising model and other clas- 
sical spin systems llZnla. In a cluster algorithm entire regions of spatially correlated 
variables are "flipped" together. This can have spectacular results near the contin- 
uum limit (or more generally a second order phase transition) in reducing critical 
slowing down. In the standard Wilson formulation of gauge theories, no such cluster 
algorithm has ever been found for models with a continuous gauge group. For the 
U(l) quantum link model, we have succeed in formulating and implementing the 
first such cluster algorithm. 

The cluster methods for the gauge model are being developed by a close analogy 
with 2-d cluster algorithms for quantum spin systems. Fat. Quantum spin systems, 
successful cluster algorithms iaizE already been invented I13'E2I and very impressive 
result obtain in simulations rWA For example efficient cluster algorithms for a 
quantum XY models (anakigous to a U(l) gauge theory) have been studied and 
we know how to construct cj improved estimators for quantum Green's functions 
such as < o''^cr~ >. By use of an explicit basis for the transfer matrix and the 
Trotter formula, one has a quantum generalization of the Fortuin-Kasteleyn map 
onto a random cluster model. The clusters for the quantum XY model are single 
closed curves in the d-fl space (see Fig We caiLprove detailed balance and basis 
independence for the quantum matrix elements t3. Also one can formulate these 
algorithms exactly on the continuous "time" interval [0, (3] making no error due to 
the usual Trotter formula approximation. 
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Figure 4: The figure shows a typical graph in Id with L = 4, and N = 2. Every white square 
represents a two spin interaction. The graph shown here splits the lattice into three clusters, 
represented by a solid line, broken line and a dotted line. 

Repiarkably an analogous procedure has been carried out for the U(l) gauge 
theory c3. The most interesting topological difference with quantum spin system 
is that now the clusters are closed surfaces (instead of closed lines) representing 
world sheets for string states propagating in the extra time direction. As in all 
efficient cluster algorithms with improved estimators, the structure of the clusters 
should signify a new set of collective co-ordinates which play an important role in 
the dynamics. For quantum link theory these collective co-ordinates evidently are 
indicative of an underlying stringy vacuum. 

Finally we should remark that in the general formulation of a particular D- 
theory, we expect that there are many alternative lattice Hamiltonians that are 
universally equivalent after dimensional reduction, but have different microscopic 
details. This flexibility at the microscopic level is (ultimately) a strength. Just 
as solvable 2-d spin systems, beginning with the classic work on the Ising model, 
have led to new theoretical insight (i.e. critical scaling) and new algorithms (i.e. 
percolation cluster methods), a particular formulation of quantum link models may 
lead to more efficient algorithms or more tractable analytic methods. The implica- 
tion is that detailed choices among possible universally equivalent options should 
go hand in hand with the development of analytical and numerical methods. There 
are many options that we are exploring that go beyond the examples I have had 
time to present. Nonetheless we feel there is a rather clear unified methodology in 
constructing a D-theory and ample reason for optimism that they can contribute 
new insight into the nature of asymptotically free field theories. 
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